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INTRODUCTION 


The  Department  of  Defense  (DOD)  has  numerous  facilities  engaged  in 
the  production  of  various  types  of  explosives  and  munitions  used  by  mil- 
itary services.  In  most  cases  the  production  of  ammunition  utilizes 
assembly  line  procedures.  Projectiles  pass  through  various  stages  of 
preparation:  filling  with  explosive,  fuzing,  marking,  and  packing. 

Hazardous  operations,  such  as  the  filling  of  the  projectile  case  with  an 
explosive  in  a powder  form  and  the  compaction  of  the  powder  by  hydraulic 
press,  are  accomplished  in  protective  cells  that  are  intended  to  confine 
the  effects  of  an  accidental  explosion. 

Most  of  the  existing  production  facilities  were  built  in  the  1940s. 
With  few  exceptions,  the  manufacturing  technology  and  existing  equipment 
represent  the  state-of-the-art  at  that  time.  The  production  equipment 
was  operated  extensively  during  World  War  II,  again  during  the  Korean 
conflict,  and  recently  during  the  Southeast  Asia  war.  Much  of  this 
equipment  and  the  housing  structures  have  been  operating  beyond  their 
designed  capacities  (Ref  1). 

DOD  is  conducting  an  ammunition  plant  modernization  program  (Ref  2) 
intended  to  greatly  enhance  safety  in  the  production  plants  by  pro- 
tective construction,  automated  processing,  and  reduction  of  the  number 
of  personnel  involved  in  hazardous  operations. 

In  1969  a joint-service  manual  (Ref  3)  was  published  to  provide 
guidance  to  the  structural  designers  of  munition  plants.  The  objectives 
of  the  manual  were  to  establish  design  procedures  and  construction 
techniques  (1)  to  prevent  propagation  of  explosions  from  one  building 
(or  part  of  a building)  to  another,  (2)  to  prevent  mass  detonations,  and 
(3)  to  protect  personnel  and  equipment.  The  manual  establishes  blast- 
load parameters  for  designing  protective  structures,  provides  methods 
for  calculating  the  dynamic  response  of  concrete  walls,  and  establishes 
construction  details  for  developing  required  strength.  The  design 
method  accounts  for  close-in  effects  of  a detonation  with  its  associated 
high  pressures  and  nonuniformity  of  loading  on  protective  barriers.  A 
detailed  method  for  assessing  the  degree  of  protection  afforded  by  a 
protective  facility  did  not  exist  prior  to  this  manual's  publication; 
consequently , the  manual  represents  a significant  improvement  in  design 
methods.  The  simplifications  made  in  the  development  of  the  design 
procedures  have  been  presented  in  the  manual.  The  analysis  of  a struc- 
ture using  the  design  procedure  will  generally  result  in  a conservative 
estimate  of  the  structure's  capacity;  therefore,  structures  designed 
using  these  procedures  will  generally  be  adequate  for  blast  loads  ex- 
ceeding the  assumed  load  conditions  (Ref  3;. 
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Ever.  with  Clio  simplifications  presented  in  Reference  3,  the  compu- 
tational procedures  are  complex  and  time-consuming . An  automated  pro- 
cedure was  required  to  give  structural  designers  the  capability  of 
performing  rapid  analysis  of  the  structural  safety  of  blast-resistant 
walls  and  doors.  The  design  parameters  interact  in  a complex  way  since 
the  procedure  is  both  nonlinear  and  dynamic.  From  a design  point  of 
view  an  optimization  procedure  was  required  to  minimize  cost  and  maxi- 
mize safety  since  blast-resistant  construction  has  been  reported  to  cost 
three  to  five  times  as  much  as  conventional  construction.  Therefore, 
the  first  objective  was  to  automate  the  analysis  procedures  for  deter- 
mining the  structural  response  of  plates  having  a bilinear  stiffness 
representation  and  subjected  to  blast  shock  and  gas  pressures.  Plates 
are  the  basic  elements  forming  sidewalls,  roofs,  floors,  and  doors  ‘■f 
cells  designed  to  confine  the  effects  of  accidental  explosions.  The 
second  objective  was  tc  provide  an  optimum  design  procedure  that  will 
automatically  produce  a least-cost  design  for  a given  geometry,  material 
properties,  and  explosive  weight  for  both  feasible  and  nonfeasible 
starting  points. 

COMPUTER  PROGRAM 

The  computer  program  was  written  in  FORTRAN  IV  for  use  with  Control 
Data  6600  series  computers.  The  program  is  composed  of  four  areas: 

1.  BJast  Load  Determination 

2.  Structural  Analysis  Parameters 

3.  Dynamic  Response 

4.  Optimization 

The  blast-] oad  determination  is  accomplished  by  subroutines  BLA, 
PIC,  SCRID,  HBA,  RATIO,  GRID,  Ga's  INTERP,  EQUIV,  HEDATA,  ARDC,  SHOCK, 
and  TNT.  The  subroutines  read  the  explosive  weight  and  type  and  cell 
geometry,  and  th»*.  compute  the  equivalent  .spherical  weight  of  TNT  and 
the  equivalent  pressure  loading  using  the  geometry  of  the  wall  and 
charge  location.  Both  the  shock  pressure  and  its  duration  and  the  gas 
pressure  and  its  duration  are  calculated  as  in  References  3 and  4. 

Using  the  duration  and  pressure  data  for  both  shock  and  gas,  the  program 
computes  an  equivalent  triangular  pressure  loading  for  each  part  and 
adds  both  together  to  produce  the  resultant  shown  in  Figure  1.  The 
total  impulse  is  then  determined  as  in  Reference  3. 

The  structural  analysis  is  accomplished  by  subroutines  SSTIFF,  DOOR 
1,  DOOR  2,  DOOR  3,  DOOR  4 and  DOOR  5.  These  routines  compute  the  stiff- 
ness, resistance,  and  equivalent  mass  of  the  plate  using  input  material 
properties  as  in  Reference  3.  Both  i lexure  and  shear  are  considered. 
Openings  in  plates  are  allowed  as  indicated  in  Figure  2c. 
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The  dynamic  response  calculation  is  accomplished  in  subroutine 
RESP.  The  program  determines  the  response  of  the  plate  modeled  as  an 
equivalent  dynamic  single-degree-of-f reedom  system  with  bilinear  stiff- 
ness and  the  pressure  loading  shown  in  Figure  1 . The  solution  technique 
is  based  on  a Newmark  iterati-  n method. 

When  a thickness  of  sand  is  specified  for  composite  construction 
(i.e.,  two  plates  with  sandfill),  the  program  computes  the  impulse 
capacity  of  the  first  plate  using  half  the  mass  of  the  sand  as  acting 
with  the  wall  as  in  Reference  3.  Figures  6-38  and  6-39  of  Reference  3 
give  the  attenuation  of  the  blast  wave  on  sand  for  evaluation  of  the 
impulse  capacity  of  the  second  wall. 

The  optimization  of  an  initial  design  is  accomplished  in  subroutines 
OPT,  MINIMZ,  PMINZ , DMXNZ,  GETE,  SUMRY , TLEFT,  and  GCOMP.  The  methodol- 
ogy used  is  that  of  a penalty  function  with  individual  minimization 
sequences  being  accomplished  by  the  Powell  method  (References  4,5,6). 


PROGRAM  INPUT 


The  program  input  consists  of  five,  or  six  cards  per  case.  Addi- 
tional cases  can  be  grouped  together.  Two  blank  cards  are  used  after 
the  last  case.  The  user's  guide,  contained  in  the  program  with  comment 
cards,  is  given  here  to  assist  in  understanding  the  input.  Card  format 
is  8F10.0  except  as  noted.  Figure  2a  is  an  input  data  sheet  to  be  used 
in  conjunction  with  Figures  2b  and  c,  which  show  the  slab  geoi.  etty  and 
orientation  that  must  be  followed.  The  input  required  for  each  card  is 


described  below. 


CARD  1 


COL 

2 

COL 

68 

COI 

69 

COL 

79 

COL 

71 

COL 

72 

COL 

73 

COI, 

74 

C.OL 

75 

COL 

76 

COL 

77 

COL 

78 

COL 

79 

COL 

80 

CARD 

2 

COL 

1 

COL 

10 

COL 

11 

COL 

20 

COL 

21 

COL 

30 

COL 

31 

COL 

40 

HEADING 

OPTIMIZATION  0 = NO  OPTIMIZATION  CALCULA- 
TION, 1 = OPTIMIZATION  CALCULATION 
FLAG  1=0  FOR  PRESSURE  CALCULATION, 

= 1 FOR  INPUT  PRESSURE  (see  Card  3) 
FLAG  2 FOR  TS  OR  Z:  0 = TS,  1 = INPUT  Z 
FLAG  3 FOR  IMPULSE  GRID:  0 = OMIT, 

1 = GRIP 

FLAG  4 0 = NO  DOOR,  1 = DOOR 

FLAG  5 PRINT:  DOOR  EQUILIBRIUM  ITERATION 

0 = OMIT,  1 = PRINT 


WEIGHT  OF  ACTUAL  EXPLOSIVE,  T.B 
EXPLOSIVE  NUMBER,  SEE  TABLE  1 
EXPLOSIVE  LENGTH /DIAMETER  RATIO 
PROJECTILE  CASE  WEICHT/EXPLOSIVE  WEIGHT 
RATIO 
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COL 

41 

COL 

50 

AMBIENT  PRESSURE  PSIA  (DEFAULT  14.69  l’SI) 

COL 

51 

COL 

60 

AMBIENT  TEMPERATURE,  °C  (DEFAULT  20°) 

COL 

61 

COL 

70 

ALTITUDE  KFT  (WHEN  PRESSURE  AND  TEMPER- 
ATURE NOT  SPECIFIED) 

COL 

71 

COL 

80 

EFFECTIVE  IMPULSE  FRACTION  COMPOSITE 
CONSTRUCTION  (see  Ref  3) 

CARD  3 

COL 

1 

COL 

10 

RA  DISTANCE  CHARGE  TO  WALL  F’T  OR  EQUAL 
IMPULSE  PSI-MS  IF  FLAG  1-1.0 

COL 

11 

COL 

20 

H WALL  HEIGHT,  FT 

COI. 

21 

COL 

30 

EL  WALL  LENGTH,  FT 

COL 

31 

COL 

40 

HLIT  HEIGHT  CHARGE  FT  OR  EQUAL  PRESSURE 
PSI  IF  FLAG  1 = 1.0 

COL 

41 

COL 

50 

EL1.IT  DISTANCE  CHARCE  TO  LEFT  SIDE 

WALL  FT  , 

COL 

51 

COL 

60 

CELL  VOLUME  FOR  GAS  PRESSURE,  Ft 

COL 

61 

COL 

70 

CELL  VENT  Ai  SA  FOR  CAS  PRESSURE,  FT 

COE 

71 

EQ  1 FOR  FLUOR  REFLECTION 

COL 

7 2 

EQ  1 FOR  Rf  OF  REFLECTION 

COL 

73 

EQ  1 FOR  I EFT  WALL  REFLECTION 

COL 

74 

EQ  1 FOR  RIGHT  WALL  REFLECTION, 
OTHERWISE,  EQ  0 FOR  NO  REFLECTION 

CARD  4 

COL 

1 

COL 

10 

DYNAMIC  YIELD  STRESS,  PSI 

COL 

11 

COL 

20 

PLATE  THICKNESS,  IN. 

COL 

21 

COL 

30 

NS  IDE  NUMBER  OF  SIDES  WALL  FIXED 

1.0  BOTTOM  SIDE  FIXED 

2.0  BOTTOM  AND  SIDE  FIXED 

3.0  2 SIDES  AND  BOTTOM  FIXED 

4.0  4 SIDES  FIXED 

5.0  SIMPLE  SUPPORTED  BEAM  AT  TOP 
AND  BOTTOM 

6.0  FIXED  BEAM  AT  TOP  AND  BOTTOM 

7.0  BEAM  BOTTOM  FIXED  TOP  SIMTLE 

13.0  3 SIDES  SIMPLE,  1 SIDE  FREE 

14.0  4 SIDES  SIMPLE 


COL 

31 

COL  40 

PLATE  HEIGHT  IF  NOT  EQUAL  TO  H CARD 

3, 

FT 

COL 

41 

COL  50 

PLATE  WIDTH  IF  NOT  EQUAL  TO  EI.  CARD 

3, 

FT 

COL 

51 

COI.  60 

ALLOWABLE  DUCTILITY  LIMIT  FOR  OPTIMIZATION 

COL 

61 

COL  70 

THICKNESS  SAND,  FT 

COL 

71 

COL  80 

E MODULUS  OF  ELASTICITY,  PSI 

4 


-IP* 


‘ITS 


>£***  norm- 


i.W**’**' 
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CARD  5 

IF  OPTION  = 1 

ON  CARD  1 

COLUMN 

73-74,  OTHERWISE  SKIP 

C01,  1 

COL  10 

Z HORIZONTAL  SECTION  MODULUS/ IN . , IN.'VlN 

COL  11 

COL  20 

7.  VERTICAL  SECTION  MODULUS/ IN.  , TN.3/IN. 

COL  21 

COL  30 

AVERAGE  MOMENT  INERTIA/IN.,  IN.4/1N. 

COL  31 

COL  40 

DOOR  WEIGHT  TOTAL,  LB 

CARD  6 

BLAST  DOOR  PARAMETERS 

IF  OPTION  = 1 ON  CARD  1 

COLUMN 

77-78,  OTHERWISE  SKIP 

COL  1 

COL  10 

DOOR  HEIGHT,  l'T 

COL  11 

COL  20 

DOOR  WIDTH,  FT 

COL  21 

COL  30 

DISTANCE  FROM  LEFT  SIDE  TO  DOOR,  FT 

COL  31 

COL  40 

DOOR  REACTION,  LB/ IN. 

OR 

COL  41 

COL  50 

DOOR  RESISTANCE  FOR  CALCULATION  OF 

COL  51 

COL  60 

REACTION,  PSI 

DISTANCE  TO  FLOOR,  FT 

NOTE:  All  values  are  fixed  point,  except  for  reflection  code  and  options 

The  explosive  number  (Card  2)  refers  Co  the  list  of  explosives  in 
Tabic  1.  This  is  used  to  compute  explosive  equivalence.  The  length./ 
diameter  ratio  for  an  explosive  sphere  is  0.0,  which  gives  a shape 
factor  of  1.0.  For  an  uncased  explosive  the  case  explosive  weight  ratio 
is  0.0.  For  sea  level  calculations,  the  ambient  air  pressure  Pai!1b,  and 
temperature  Taln^,  and  altitude  can  be  left  blank  and  will  default  to 
14.69  psi  and  20°C.  If  the  flag  in  the  heading  card  is  set  to  1,  the 
impulse,  duration,  and  pressure  will  be  read  on  Card  3.  If  the  flag  is 
left  blank,  the  charge  to  wall  distance,  charge  height,  and  distance 
from  the  left  side  will  be  read.  If  NSIDE  is  left  blank,  the  program 
will  sum  the  number  of  reflecting  sidewall  surfaces  specified  on  Card  3. 
The  separate  use  of  NS ll»F  is  helpful  when  a frangible  wall  is  present, 
which  creates  a shock  reflection  but  does  not  provide  any  support. 

When  optimization  and  composite  construction  are  specified  together, 
the  program  will  optimize  the  design  to  resist  the  given  or  computed 
impulse.  For  the  cas,  when  two  walls  are  acting  together-each  resisting 
a portion  of  the  impulse— it  is  necessary  to  specify  the  effective  impulse 
to  be  applied  to  the  wall  under  design.  The  total  impulse  is  multiplied 
by  the  decimal  number  specified  on  Card  2.  This  procedure  is  based  on 
similar  work  for  concrete  (Ref  4). 

The  NS1DF.  (see  Figure  2b)  conditions  1 through  4 are  intended  to  be 
used  to  represent  steel  cell  walls  and  roofs;  NS IDE  conditions  3 through 
7 are  steel  plates  spanning  in  one  direction.  The  NS IDE  conditions  13 
and  14  are  specifically  intended  to  represent  typical  steel,  plate  doors 
and  pass-through  windovis. 
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The  optimization  problem  consists  of  finding  the  least-cost  struc- 
ture that  satisfies  all  the  design,  constraints;  or,  stated  in  opti- 
mization terms:  Find  X such  that  M(X)  is  a minimum  and 

g.(X)  <0  i = 1,  2,  N 

-V 

where  X = vector  of  design  variables 

N = number  of  design  constraints 

g = vector  of  design  constraints 

M = objective  function 

Specifically  for  this  problem,  the  design  variables  selected  are  areas 
of  steel  reinforcement  and  thickness  of  concrete.  The  design  constraints 
are  the  flexural  and  shear  limits.  The  objective  function  consists  of 
the  costs  of  formwork  and  concrete  flexural  and  shear  reinforcement. 

Fixed  Variables 

W = explosive  weight 

H = height 
HI.  = length 

h = height  of  explosive  above  floor 

i = distance  of  explosive  from  left  side  of  wall 

= distance  of  explosive  from  wall 
1 = reflection  code 

i = dynamic  yield  stress 

P = ductility 

Design  Parameters,  X 

X = t (thickness  of  plate) 

Constraints,  g(X) 

6 (X)  = 6(0)  , maximum  deflection 

t > 0.05  minimum  thickness 

t < 20  maximum  thickness 
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Tin.’  methodology  (Kef  i>  and  6)  selected  uses  Die  unconstrai  nod 
in  i uimizat i on  approach.  The  problem  is  converted  to  an  unconstrained 
minimization  by  construct  inn  a function  4s  of  the  general  form 

4(X,r)  = M(X)  + HsijOb,  Kn(X),  rj 

For  this  problem  tlie  Interior  penalty  function  technique  was  selected. 
This  methodology  is  suitable  when  gradients  are  not  available,  and, 
because  the  method  uses  the  feasible  region,  a usable  solution  always 
results.  The  objective  function  is  augmented  with  a penalty  term  that 
is  small  at  points  away  from  the  constraints  in  the  feasible  region  but 
increases  rapidly  as  the  constraints  are  approached.  The  form  is  as 
f ol  lows: 


. 

®(X.r)  = M(X)  - r £ - 

j*l  Rj(X) 

where  M is  to  be  minimized  over  all  X satisfying  g (X)  ^ 0,  j = 2 ...N. 
Note  that  if  r is  positive,  then,  since  at  any  interior  point  all  ol  the 
terms  in  the  sum  arc  negative,  the  effect  is  to  add  a positive  penalty 
to  M(X).  As  the  boundary  is  approached,  some  g (X)  will  approach  zero, 
and  the  penalty  will  increase  rapidly.  The  parameter,  r,  will  he  made 
successively  smaller  in  order  to  obtain  the  constrained  minimum  of  M. 

Objective  unction,  F 

Cost  = F = H • EL  • t • C 

where  C = volumetric  cost  of  material 


where  r = penalty  parameter. 

The  program  requires  a starting  point  in  the  feasible  region  before 
optimization  can  proceed.  This  is  accomplished  automatically  by  the 
program  by  incrementing  the  design  variables  until  a feasible  point  is 
reached . 

An  algorithm  which  comprises  the  steps  most  commonly  used  is  as 
f ol lows : 


1.  Given  a starting  point  XQ,  satisfying  all  g^(X)  < 0,  and  an 
initial  value  for  r,  minimize  to  obtain  XD  . 

2.  Check  for  convergence  of  X . to  the  optimum. 

mm 


7 


3.  If  the  convergence  criterion  is  not  satisfied,  reduce  r by 
r ■*-  rc,  where  c < 1. 

4.  Compute  a new  starting  point  for  the  minimization,  initialize 
the  minimization  algorithm,  and  repeat  from  step  1. 

The  logic  diagram  for  the  interior  penalty  functions  technique  is 
shown  in  Figure  3. 

The  minimization  for  4>(X,r)  shown  in  Figure  3 is  accomplished  by  a 
method  developed  by  Powell  using  conjugate  directions  (Ref  5 and  6). 

Powell's  method  can  be  understood  as  follows:  Given  that  the 
function  has  been  minimized  once  in  each  of  the  coordinate  directions 
and  then  in  the  associated  pattern  direction,  discard  one  of  the  coor- 
dinate directions  in  favor  of  the  pattern  direction  for  inclusion  in  the 
next  m minimizations,  since  this  is  likely  to  be  a better  direction  than 
the  discarded  coordinate  direction.  After  the  next  cycle  of  minimizations, 
generate  a new  pattern  direction  and  again  replace  one  of  the  coordinate 
directions.  This  process  is  illustrated  in  Figure  4. 

Figure  5 is  a logic  diagram  for  the  unconstrained  minimization 
algorithm.  The  pattern  move  is  constructed  in  block  A,  then  used  for  a 
minimization  step  (blocks  B and  C) , and  then  stored  in  Sn  (block  D)  as 
all  of  the  directions  are  up-numbered  and  S]_  is  discarded.  The  direc- 
tion Sn  will  then  be  used  for  a minimizing  step  just  before  the  construc- 
tion of  the  next  pattern  direction.  Consequently,  in  the  second  cycle, 
both  X and  Y in  block  A are  points  that  are  minima  along  Sn,  the  last 
pattern  direction.  This  sequence  will  impart  special  properties  to  Sn+i 
= X - Y that  are  the  source  of  the  rapid  convergence  of  the  method. 

Figure  5 shows  a^block  requiring  a one-dimensional  minimization  of 
a*  of  the  function  (f)(X  + a Sq).  The  one-dimensional  minimization  uses  a 
four-point  cubic  interpolation.  It  finds  the  minimum  along  the  direc- 
tion Sq,  where  X is  the  coordinate  of  the  previous  minimum.  By  trial 
and  error  it  finds  three  points  with  the  middle  one  less  than  the  other 
two.  It  makes  a quadratic  interpolation  and,  then,  a cubic  interpolation. 
If  the  actual  function  evaluated  at  the  new  interpolated  point  is  not 
sufficiently  close  to  that  of  the  preceding  point,  or  if  it  is  not 
sufficiently  close  to  the  interpolated  function,  then  another  cubic 
interpolation  is  made.  The  logic  for  this  algorithm  is  shown  in  Figure 
6. 


APPROXIMATE  COMPUTATION  OF  DOOR  REACTION 

It  should  be  emphasized  that  this  program  is  intended  to  assist  in 
rapid  approximate  design  and  not  detailed  analysis.  The  basic  procedures 
in  References  3 and  4 and  used  herein  have  been  found  to  be  sufficiently 
accurate  for  simple  geometries  of  beams  and  slabs  without  openings. 

Figure  7 compares  deflections  for  a plate  fixed  on  four  sides  and  for  a 
beam;  the  approximate  solutions  and  the  finite  element  solutions  agree 


8 


within  about  10%.  However,  the  static  shear  procedures  suggested  in 
Reference  3 are  seen  to  be  substantially  below  dynamic  shears  (figure 
8);  this  is  a limitation  of  the  approximate  procedures  and  is  under 
current  investigation. 

A steel  door  attached  to  a concrete  wall  was  examined  using  a 
finite  element  technique.  Figure  9 shows  the  slab  and  door;  Figure  10 
shows  the  deflection  of  the  door  by  the  approximate  procedure  developed 
herein  and  the  finite  element  procedure.  There  is  some  disagreement  in 
deflection,  especially  when  one  considers  the  deflecting  top  support. 

It  should  be  particularly  noted  that  the  deflecting  support  condition 
for  actual  doors  on  slabs  (modeled  correctly  by  finite  element  and 
assumed  rigid  by  approximate  solution)  absorbs  significant  amounts  of 
energy  by  rigid-body/door  motion.  Thus,  the  resulting  center  door  de- 
flection is  reduced.  The  resulting  dynamic  shear  around  the  door  (trans- 
ferred to  the  wall)  is  reduced  from  what  would  be  computed  for  a nonde- 
flecting plate  using  approximate  dynamic  plate  theory  (Figure  11).  The 
alternatives  are  to  use  finite  element  analysis  procedures  or  to  modify 
dynamic  plate  theory.  Finite  element  analysis  is  certainly  the  better 
approach;  however,  it  is  basically  an  analysis  technique  and  is  more 
difficult  and  expensive  to  use  than  the  simpler  approximate  procedure. 

Tt  is  suggested  that  the  shear  calculated  from  approximate  plate  theory 
be  adjusted  by  a constant  for  use  as  a door  reaction  required  for  input 
to  wall  design  (Ref  4). 

The  maximum  reaction  occurs  at  the  moment  the  slab  first  reaches 
yield.  At  this  point  the  combination  of  load  and  resistance  is  maximum. 
Table  'l  gives  maximum  dynamic  reaction  for  a simply  supported  plate. 

For  the  case  of  one  side  free  and  three  sides  simple-supported,  the  b- 
dimension  doubled  may  be  used.  The  values  of  pressure  P and  resistance 
R are  taken  from  the  computer  output  3t  the  time  of  yielding.  The 
reaction  values  should  be  adjusted  for  support  deflection.  The  value  of 
1.0  is  suggested  for  nondeflecting  supports  and  0.5  for  full  deflecting 
supports  as  approximate  factors.  Once  design  has  been  finalized  it  is 
suggested  that  results  be  analyzed  using  a finite  element  analysis. 


DISCUSSION 


This  program  was  developed  to  perform  rapid  design  of  steel  plates 
used  to  form  the  sides  and  roofs  of  blast  cells  and  also  of  steel  plates 
used  as  doors.  Provisions  are  included  for  use  of  plastic  section 
modules  for  built-up  doors;  but  optimization  of  such  doors  may  not  be 
performed  because  the  weight-strength  function  is  not  defined. 

In  general,  the  methods  used  in  th“  coi  *uter  program  follow  Reference 
3;  consequently,  the  accuracy  of  both  in  ’.nr  same.  These  are  discussed 
in  detail  in  References  3 and  4 and  will  not  be  presented  here.  The 
solution  of  the  dynamic  response  equation  of  motion  has  been  found  to 
agree  very  closely  with  the  response  chat t of  Reference  3.  Additionally, 
the  solution  covers  a wider  range  and,  thus,  is  more  accurate  in  the 
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areas  not  defined  by  the  response  chart.  When  the  loading  is  less  than 
one  hundredth  of  the  natural  period,  the  response  is  determined  by 
impulse  equilibrium.  The  basic  dynamic  model  is  limited  to  the  primary 
response  mode  and  does  not  consider  higher  modes. 

The  blast  impulse  computation  is  restricted  to  a geometry  in  which 
the  slab's  height-to-length  ratio  is  greater  than  0.2.  A modification 
was  made  by  the  Naval  Surface  Weapons  Center  to  the  original  Picatinny 
Arsenal  Program  to  remove  several  minor  problem  areas,  such  as  the 
location  of  the  charge.  The  blast  impulse  has  all  the  limitations 
associated  with  the  original  Picatinny  programs  that  are  caused  by 
limitations  in  the  test  data.  It  assumes  the  charge  is  an  equivalent 
sphere  of  TNT.  Shape  effects,  explosive  equivalence,  and  explosive 
casings  are  considered,  but  only  in  an  empirical  manner  as  a result  of 
limited  available  data. 

Example  problems  are  presented  in  the  Appendix. 
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Table  1.  List  of  Explosives 


Explosive 

Number 


Explosive  Name  and  Composition 


1 


3 

4 

5 

6 

7 

8 
9 

10 


11 


12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 
29 


TNT 

TNETB 

EXPLOSIVE  D 

PENTOLITE  (PETN/TNT  50/50) 
PICRATOL  (EXPLOSIVE  D/TNT  52/48) 
CYCLOTOL  (RDX/TNT  70/30) 

COMP  B (RDX/TNT/WAX  59.4/39.6/1.0) 
RDX/WAX  (98/2) 

COMP  A-3  (RDX/WAX  91/9) 
TNETB/AL  (90/10) 

TNETB/AL  (78/22) 

TNETB/AL  (72/28) 

TNETB/AL  (65/34) 

TRITONAL  (TNT/AL  80/70) 
RDX/AL/WAX  (88/10/2) 
RDX/AL/WAX  (89/20/2) 
RDX/AL/WAX  (74/21/5) 
RDX/AL/WAX  (74/22/4) 
RDX/AL/WAX  (62/33/5) 

TORPEX  II  (RDX/TNT/AL  42/40/18) 

H6  (RDX/TNT/AL/WAX  45/29/21/5) 
HBX-1  (RDX/TNT/AL/WAX  40/38/16/5) 
HBX-3  (RDX/TNT/AL/WAX  31/29/35/5) 
TNETB /RDX/AL  (39/26/35) 
ALUMINUM 
WAX 
RDX 
PETN 
TETRYL 
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Table  2.  Four  Sides,  Uniform  Load* 


Simple 

support 


Strain 

a/b 

Dynamic  Reactions** 

Range 

V/b 

A 

Va 

1.0 

0.07P  + 0.18R 

0.07P  + 0.18R 

0.9 

0.06P  + 0.16R 

0.08P  + 0.20R 

Elastic 

o 

CO 

0.06P  + 0.UR 

0.08P  + 0.22R 

0.7 

0.05P  + 0.13R 

C.08P  + 0.24R 

0.  6 

0.04P  + 0.11R 

0.09P  + 0.26R 

0.  5 

0.04P  + 0.09R 

0.09P  + 0.28R 

1.0 

0.09P  + O.lbR 

m 

0. 09 P 4-  0. 16K 

m 

0.9 

0.08P  + 0.15R 

m 

0.09P  + 0.18R 

m 

Plastic 

O 

00 

0.07P  + 0.13R 

m 

0.10P  + 0.20R 

m 

0.7 

0.06P  + 0.12R 
m 

0.10P  + 0. 22R 

m 

0.6 

0.05P  + 0.10R 

m 

0.10P  + 0 . 25R 

ip 

0.5 

0. GAP  + 0.08R 
m 

0.11P  + 0.27R 

TP 

*Rased  on  information  from  Ref  7. 


**P  = pressure  at  time  of  yield,  psi 

R = resistance,  psi 
Rm  = yield  resistance,  psi 


Start 


Figure  5.  Logic  diagram  for  minimization  of  4>(X) 


CEL 


exit 


satisfies  I'j  > l;,  > l'2  ”r  1 1 F3  > V2 
Figure  6.  One-dimensional  minimization  algorithm. 
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level  of  shear  indicated  Ret  3 


1 Ib/in.  = 0.179  kg/cm 


(in.) 


0.9 


approximate  <T. 


3.0  rns) 


24 


APPENDIX 


EXAMPLE  PROBLEM  1 

Design  a door  made  of  steel  plate  f;  r the  following: 

1.  Door  Height  = 6 ft 

2.  Door  Width  = 4 ft 

3.  Dynamic  Yield  Stress  = 48,000  psi 

4.  Simple  Support,  Bottom  Free 

5.  Allowable  Ductility  = 10 

The  door  is  contained  in  a wall  12  ft  wide  by  10  ft  high.  Side  walls 
and  roof  are  present  to  provide  reflecting  surfaces.  The  explosive  is 
10  lb  Composition  B encased  located  3 ft  from  the  wall,  5 ft  from  the 
left  side,  and  3 ft  above  the  floor.  Figure  A- 1 shows  the  example 
problem  input  form,  and  Figure  A-2  shows  the  output. 

EXAMPLE  PROBLEM  2 

Design,  a steel,  plate  window: 

1.  Height  = 3 ft 

2.  Width  = 3 ft 

3.  Dynamic.  Yield  Stress  = 48,000  psi 

4.  Simple  Support,  Four  Sides 

5.  Allowable  Ductility  = 10 

The  window  is  located  on  the  wall  of  a cell  12  ft  long  by  10  ft  high. 

A 12-lb  TNT  explosive  with  length- to-diameter  of  2.5  and  case-to- 
explosivt  of  1.2  is  located  3 ft  away  from  the  wall,  5 ft  from  the  left 
side,  and  3 ft  above  ground.  The  cell  has  two  sidewalls  and  a floor. 

No  roof  reflecting  surface  is  present.  Figure  A-3  shows  the  example 
piroblem  input  and  Figure  A-4  gives  the  output. 


COHP  P (RQX/TNT/.vAX,  59, 4/39. o/  1 ,0) 

EXPLOSIVE  PROPERTIES CHARGE  *E I ?HT (LB)  * 10.00 
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Figure  A-2.  Continued 
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Figure  A-4.  Continued 
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Figure  A-4.  Continued 
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Figure  A- 'i . Continued 
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Figure  A- A.  Continued 
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Figure  A-A.  Continued 
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ARMS  PI  AST  EC'  TVatmny  Arsenal  ( A M Anralnne.  SMLPA-FR-M  D)  Dover  NJ 

ASST  SECRETARY  OF  THE.  NAVY  Spec  Assist  F.nenty  tP  Waterman).  Washington  DC 

BUREAU  OF  RFClAMATTONCodc  t<12(C  ScUnder)  Denver  CO 

(TNCI.ANT  Civil  E.ngr  Supp  Plans  Ofr  Norfolk.  S' A 

CINCPAC  F a.  Fngrng  Dis  |J44)  Makalapi.  HI 

CNM  NM  AT  OKT246 ( Dtetcrle)  Wash.  DC 

CNOCode  NOP-dW.  Washington  IXYOPMfTj  (J  Boomtuh).  Pentagon 
CT)M(X  FIANSYSPAC  SC  F.  Pearl  HarKsr  HI 

DEFENSE  CIVIL  PREPAREDNESS  AGENCY  J O.  Buchanan.  Washington  DC 
DEFENSE  DOCUMENTATION  CTR  Alexandria.  VA 
DEFENSE  INTELLIGENCE  AGENCY  Du  . Washington  IX 
t>NA  STTL.  W axhmgton  DC 

DOD  Explosives  Safely  Roan]  (Library ).  Washington  DC 
IX)H  Dr.  C ohen 

DTNSRIX’Code  i70h.  Bethcsda  MD;  Code  172  (St  Krcnzke).  Bethexda  MD 
DTNSRIX'  Code  4121  (R.  Riven).  Annapolis.  StD 
FT.TCOMBATTRACENl.ANT  PWO.  Virginia  Bch  VA 

MARINE  CORPS  BASE  Camp  Pendleton  CA  92055:  Code  43-260.  Camp  Lcjcunc  NC:  M A R Division,  Camp  Lejcune 
NC;  PW'O.  CampS.  D.  Butler,  Kawasaki  Japan 
MARINE  CORPS  HQS  Code  LFF-2.  W ashington  DC 

MCAS  Facil.  F.ngr.  Div.  Cherry  Point  NC;  Code  PWE.  Kaneohe  Bay  HI;  Code  S4.  Quantico  VA;  J.  Taylor.  Iwakuni 
Japan;  PWD.  Dir.  Maint.  Control  Div..  Iwakuni  Japan;  PWO  Kaneohe  Bay  HI 
MCLSBB520.  Barstow  CA 
MCRD  PWO,  San  IXcgo  Ca 

NADCodeOl  IB-1.  Hawthorne  NV;  Engr.  Dir.  Hawthorne.  NV 
NAF  PW’O  Sigonclla  Sicily:  PW’O.  Aisugi  Japan 

NAS  CO.  Guantanamo  Bay  Cuba:  Code  1 14,  Alameda  CA;  Code  183  (Fac.  Plan  BR  MGR);  Code  18700.  Brunswick 
ME:  Code  6234  (G.  Trask).  Point  MuguCA:  Dir.  Util.  Div.,  Bermuda:  ENS  Buchholz.  Pensacola,  FL:  PW  (J. 
Maguire).  Corpus  Chrisii  TX;  PWD  Maint.  Div.  . New  Orleans.  Belle  Chasse  LA:  PW0(M.  Elliott).  Los  Alamitos 
CA;  PWO  Belle  Chasse.  LA:  PW'O  Chase  Field  Beeville.  TX;  PWO  Key  West  FL;  PWO  Whiting  Fid.  Milton  FL: 
PW'O.  Dallas  TX;  PWO.  Glenview  IL;  PWO,  Kingsville  TX:  PWO.  Mir.  mar.  San  Diego  CA;  SCE  Lant  Fleet 
Norfolk.  VA:  SCE  Norfolk.  VA:  SCE.  Barbers  Point  HI 
NATL  RESEARCH  COUNCIL  Naval  Studies  Board.  Washington  DC 
NAVACT  PWO.  London  UK 
NAVAEROSPREGMEDCEN  SCE.  Pensacola  FL 
NAVAL  FACILITY  PWO.  Barbados;  PWO.  Brawdy  Wales  UK 

NAVCOASTSYSLAB  CO,  Panama  City  FL;  Code  715  (J.  Quirk)  Panama  City.  FL:  Library  Panama  City.  FL 
NAVCOMM  AREAMSTRSTA  PWO.  Norfolk  VA;  PWO.  Wahiawa  HI:  SCE  Unit  I Naples  Italy 
NAVCOMMSTACode  401  Nea  Makri. Greece;  PWO.  Adak  AK;  PWO,  Exmouth.  Australia  1 
NAVCONSTRACEN  CO  (CDR  C.L.  Neugcnt).  Port  Hueneme.  CA 
NAVEDTRAPRODEVCEN  Tech.  Library 
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NAVEODFAC  Code  605.  Indian  Head  MD 

NAVFACF.NGCOM  Code  043  Alexandria.  V A:  Code  044  Alexandria.  VA:  Code  0451  Alexandria.  V A:  Code  0453  (D. 
Polter)  Alexandria.  VA:  Code  0454B  Alexandria.  Va:  Code 04 B5  Alexandria.  V A:  Code  1023  (T.  Stevens) 
Alexandria.  VA;  Code  104  Alexandria.  VA:  Code  2014  (Mr.  Taam).  Pearl  Harbor  HI;  Morrison  Yap.  Caroline  Is.: 
PL-2  Ponce  P.R.  Alexandria.  VA 

NAVFACENGCOM  - CHES  DIV.  Cede  101  Wash.  DC;  Code  402  (R.  Morony)  Wash.  DC:  Code  FPO-I  (C.  Bodey) 
Wash.  DC:  Code  FPO-1  (Ottsen)  Wash.  DC:  Code  FPO-1 SP  (Dr.  Lewis)  Wash.  DC:  Code  FPO- 1 SPI 3 (T  F 
Sullivan)  Wash.  DC:  Code  FPO-IP12  (Mr.  Scola).  Washington  DC:  Scheessele.  Code  402.  Wash.  DC 
NAVFACENGCOM  - LANT  DIV.:  Eur.  BR  Deputy  Dir.  Naples  Italy:  RDT&ELO  09P2.  Norfolk  VA 
NAVFACENGCOM  - NORTH  DIV.  Code  1028.  RDT&ELO.  Philadelphia  PA:  Design  Div.  (R.  Masino).  Philadelphia 
PA:  ROICC.  Contracts.  Crane  IN 

NAVFACENGCOM  - PAC  DIV.  Code  402.  RDT&E.  Pearl  Harbor  HI:  Commander.  Pearl  Harbor.  HI 
NAVFACENGCOM  - SOUTH  DIV.  Code  90.  RDT&ELO.  Charleston  SC:  Dir..  New  Orleans  LA 
NAVFACENGCOM -WEST  DIV.  Code04B:  O9P/20:  RDT&ELO  Code  201 1 San  Bruno.  CA 
NAVFACENGCOM  CONTRACT  AROICC.  Point  Mugu  CA:  AROICC.  Quantico.  VA:  Eng  Div  dir.  Southwest  Pac. 
Manila.  Pi:  OICC.  Southwest  Pac.  Manila.  PI:  OICC/ROICC.  Balboa  Canal  Zone:  ROICC  (Ervin)  Puget  Sound 
Naval  Shipyard.  Bremerton.  WA;  ROICC  AF  Guam:  ROICC  LANT  DIV..  Norfolk  VA:  ROICC  Off  Point  Mugu. 
CA:  ROICC.  Keflavik.  Iceland:  ROICC.  Pacific.  San  Bruno  CA 
NAVMAG  SCE.Guam 
NAVMIRO  OIC.  Philadelphia  PA 

NAVNUPWRU  MUSE  DET  Code  NPU80(ENS  W.  Morrison).  Port  Hueneme  CA 
NAVOCEANSYSCEN  Code  6565  (Tech.  Lib.).  San  Diego  CA;  Research  Lib..  San  Diego  CA 
NAVORDSTA  PWO.  Louisville  KY 
NAVPETOFF  Code  30.  Alexandria  VA 
NAVPGSCOL  Code  61 WL  (O.  Wilson)  Monterey  CA 

NAVPHIBASE  CO.  ACB  2 Norfolk,  VA;  Code  S3T.  Norfolk  VA:  Harbor  Clearance  Unit  Two.  Little  Creek.  VA 
NAVREGMEDCEN  Code  3041 . Memphis.  Millington TN:  SCE  (D.  Kaye):  SCE (LCDR  B.  E.  Thurston).  San  Diego 
CA:  SCE.  Camp  Pendleton  CA 

NAVSCOLCECOFF  C35  Port  Hueneme.  CA:  C44A  (R.  Chittenden),  Port  Hueneme  CA:  CO.  Code  C44A  Port 
Hueneme.  CA 

NAVSEAS YSCOM  Code  OOC  (LT  R.  MacDougal).  Washington  DC 
NAVSEC  Code  6034  (Library).  Washington  DC 
NAVSECGRU  ACT  PWO,  Torn  Sta,  Okinawa 
NAVSHIPREPFAC  Library.  Guam;  SCE  Subic  Bay 

NAVSHIPYD;  Code  202.4.  Long  Beach  CA:  Code  202.5  (Library)  Puget  Sound.  Bremerton  WA:  Code  400,  Puget 
Sound:  Code  404  (LT  J.  Riccio),  Norfolk,  Portsmouth  VA;  Code  410.  Mare  Is..  Vallejo  CA:  Code  440  Portsmouth 
NH;  Code  440.  Norfolk;  Code  440,  Puget  Sound.  Bremerton  WA:  Code  440.4.  Charleston  SC;  Library.  Portsmouth 
NH;  PWD(LTN.B.  Hall).  Long  Beach  CA:  PWO.  Mare  Is.;  Tech  Library.  Vallejo.  CA 
NAVSTA  CO  Naval  Station.  Mayport  FL;  CO  Roosevelt  Reads  P.R.  Puerto  Rico:  Engr.  Dir..  Rota  Spain:  Maim.  Div. 
Dir/Code  531,  Rodman  Canal  Zone:  PWD(LTJG.P.M.  Motolenich),  Puerto  Rico:  PWO.  Keflavik  Iceland:  PWO. 
Mayport  FL;  SCE,  Guam:  SCE,  Subic  Bay.  R.P.;  Utilities  Engr  Off.  (LTjG  A.S.  Ritchie).  Rota  Spain 
NAVSUBASE  LTJG  D.W.  Peck.  Groton.  CT 

NAVSUPPACT  CO.  Seattle  WA;  Code  413.  Seattle  WA;  Engr.  Div.  (F.  Mollica).  Naples  Italy:  LTJG  McGarrah. 
Vallejo  CA 

NAVSURFWPNCEN  PWO.  White  Oak.  Silver  Spring.  MD 
NAVTECHTRACEN  SCE.  Pensacola  FL 

NAVWPNCEN  Code  2636  (W.  Bonner),  China  Lake  CA;  PWO  (Code  26).  China  Lake  CA:  ROICC  (Code  702).  China 
Lake  CA 

NAVWPNSTA  ENS  G.A.  Lowry,  FallbrookCA:  Maint.  Control  Dir..  Yorktown  VA:  PW  Office  (Code09Cl) 
Yorktown,  V A:  Security  Offr,  Earl.  Colts  Neck  NJ 
NAVWPNSUPPCEN  Code  09  (Boennighausen)  Crane  IN 
NCBU  405  OIC.  San  Diego,  CA 

NCBC  CEL  (CAPTN.  W.  Petersen),  Port  Hueneme,  CA:  CEL  A01C  Port  Hueneme  CA;  Code  10  Davisville,  RI; 

Code  155,  Port  Hueneme  CA:  Code  156.  Port  Hueneme.  CA:  PW  Engrg.  Gulfport  MS:  PWO  (Code  80)  Port 
Hueneme,  CA 
NCR  20,  Commander 

NMCB  133  (ENST.W.  Nielsen):  5.  Operations  Dept.;  74,  CO;  Forty,  CO:  THREE,  Operations  Off. 

NORDA  Code  440  (Ocean  Rsch.  Off)  Bay  St.  Louis,  MS 

NRLCode  8400  (J.  Walsh),  Washington  DC:  Code  8441  (R.A.  Skop),  Washington  DC 
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NSC  Code  54.1  (Wynne),  Norfolk  VA 
NSDSCE.  Subic  Bay.  R.P. 

NTCCcide  54  (ENS  P.  G.  Jacket).  Orlando  FL:  Commander  Orlando.  FL 

NUSCCode  131  New  London.  CT;  Code  EA123  (R.S.  Munn).  New  London  CT:  Code  TA 131  (G.  De  la  Cruz).  New 
London  CT 

ONR  Code  700F  Arlington  VA;  Dr.  A.  Laufer,  Pasadena  CA 
PMTC  Code  4253-3.  Point  Mugu.  CA;  Pat.  Counsel.  Point  Mugu  CA 

PWC  ENS  J.E.  Surash,  Pearl  Harbor  HI;  ACE  Office  (LTJG  St.  Germain)  Norfolk  VA:  CO  Norfolk.  VA:  CO.  Great 
Lakes  IL;  Cede  1 16  (LTJG.  A.  Eckhart)  Great  Lakes.  IL;  Code  120.  Oakland  CA:  Code  1 20C( Library)  San  Diego. 
CA:  Code  128.  Guam:  Code  200.  Great  Lakes  IL:  Code  200.  Oakland  CA;  Code  220  Oakland.  CA:  Code  220.1. 
Norfolk  VA:  Code  30C  (Boettcher)  San  Diego.  CA:  Code  680.  San  Diego  CA:  OIC  CBU-405.  San  Diego  CA:  XO 
Oakland.  CA 

U.S.  MERCHANT  MARINE  ACADEMY  Kings  Point.  NY  (Reprint  Custodian) 

USAF  SCHOOL  OF  AEROSPACE  MEDICINE  Hyperbaric  Medicine  Div.  Brooks  AFB.  TX 

USCG  (G-ECV)  Washington  Dc;  (G-ECV/61 ) (Burkhart)  Washington.  DC:  G-EOE-4/61  (T.  Dowd).  Washington  DC 

USCG  ACADEMY  LTN.  Stramandi,  New  London  CT 

USCG  R&D  CENTER  Tech.  Dir.  Groton.  CT 

USNA  Ch.  Mech.  Engr.  Dept  Annapolis  MD:  PWD  Engr.  Div.  (C.  Bradford)  Annapolis  MD;  PWO  Annapolis  MD 
AMERICAN  CONCRETE  INSTITUTE  Detroit  MI  (Library) 

CALIFORNIA  STATE  UNIVERSITY  LONG  BEACH.  CA(CHELAPATI) 

CORNELL  UNIVERSITY  Ithaca  NY  (Serials  Dept.  Engr  Lib.) 

DAMES  & MOORE  LIBRARY  LOS  ANGELES.  CA 

FLORIDA  ATLANTIC  UNIVERSITY  Boca  Raton  FL  (Ocean  Engr  Dept..  C.  Lin) 

FLORIDA  ATLANTIC  UNIVERSITY  Boca  Raton  FL  (W.  Tessin) 

FLORIDA  f ECHNOLOGICAL  UNIVERSITY  ORLANDO.  FL  (HARTMAN) 

GEORGIA  INSTITUTE  OFTECHNOLOGY  Atlanta  GA  (School  r* Civil  Engr..  Kahn) 

LEHIGH  UNIVERSITY  Bethlehem  PA  (Fritz  Engr.  Lab  No.  13,  L edle);  Bethlehem  PA  (Linderman  Lib.  No.30. 
Flecksteiner) 

I JBRARY  OF  CONGRESS  WASHINGTON.  DC  (SCIENCES  & TECH  DIV) 

MICHIGAN  TECHNOLOGICAL  UNIVERSITY  Houghton.  MI  (Haas) 

MIT  Cambridge  MA:  Cambridge  MA  (Rm  10-500.  Tech.  Reports.  Engr.  Lib.):  Cambridge  MA  (Whitman) 

NY  CITY  COMMUNITY  COLLEGE  BROOKLYN.  NY  (LIBRARY) 

UNIV.  NOTRE  DAME  Katona.  Notre  Dame.  IN 

OREGON  STATE  UNIVERSITY  CORVALLIS.  OR  (CE  DEPT.  HICKS) 

PENNSYLVANIA  STATE  UNIVERSITY  UNIVERSITY  PARK.  PA(GOTOLSKI) 

PURDUE  UNIVERSITY  Lafayette  IN  (Leonards);  Lafayette.  IN  (CE  Engr.  Lib) 

SAN  DIEGO  STATE  UNIV.  Dr.  Krishnamoorthy,  San  Diego  CA 
SEATTLE  U Prof  Schwaegler  Seattle  WA 

SOUTHWEST  RSCH  INSTJ.  Maison,  San  Antonio  TX;  R.  DeHart.  San  Antonio  TX 
STANFORD  UNIVERSITY  Engr  Lib,  Stanford  CA:  Stanford  CA  (Gene) 

STATE  UNIV.  OF  NEW  YORK  Buffalo.  NY 

TEXAS  A&M  UNIVERSITY  COLLEGE  STATION.  TX  (CE  DEPT):  College  Station  TX  (CE  Dept.  Herbich) 
UNIVERSITY  OF  CALIFORNIA  BERKELEY,  CA  (CE  DEPT,  GERWICK):  BERKELEY.  CA  (OFF.  BUS.  AND 
FINANCE.  SAUNDERS);  Berkeley  CA  (B.  Bresler);  Berkeley  CA  (Dept  of  Naval  Arch.);  Berkeley  CA  (R. 
Williamson);  DAVIS,  CA  (CE  DEPT.  TAYLOR);  LIVERMORE.  CA  (LAWRENCE  LIVERMORE  LAB. 
TOKARZ) 

UNIVERSITY  OF  DELAWARE  Newark.  DE  (Dept  of  Civil  Engineering.  Chesson) 

UNIVERSITY  OF  HAWAII  Dr  Chiu  Honolulu,  HI;  HONOLULU,  HI  (SCIENCE  ANDTECH.  DIV.):  Honolulu  HI 
(Dr.  Szilard) 

UNIVERSITY  OF  ILLINOIS  Metz  Ref  Rm,  Urbana  IL;  URBANA.  IL  (LIBRARY);  URBANA.  IL  (NEWARK): 
Urbana  IL  (CE  Dept.  W.  Gamble) 

UNIVERSITY  OF  MASSACHUSETTS  (Heronemus).  Amherst  MA  CE  Dept 
UNIVERSITY  OF  MICHIGAN  Ann  Arbor  MI  (Richart) 

UNIVERSITY  OF  NEBRASKA-LINCOLN  Lincoln.  NE (Ross  Ice  Shelf  Proj.) 

UNIVERSITY  OF  NEW  MEXICO  J Nielson-Engr  Matls  & Civil  Svs  Div.  Albuquerque  NM 
UNIVERSITY  OF  TEXAS  Inst.  Marine  Sci  (Library),  Port  Arkansas  TX 
UNIVERSITY  OFTEXAS  AT  AUSTIN  AUSTIN.  TX  (THOMPSON);  Austin.  TX  (Breen) 

UNIVERSITY  OF  WASHINGTON  Dept  of  Civil  Engr  (Dr.  Mattock),  Seattle  WA;  SEATTLE.  WA  (MERCHANT): 
SEATTLE.  WA  (OCEAN  ENG  RSCH  LAB.  GRAY) 
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URS  RESEARCH  CO.  LIBRARY  SAN  MATEO.  CA 

ALFRED  A.  YEE& ASSOC.  Honolulu  HI 

AMETEK  Offshore  Res.  & Engr  Div 

APPLIED  TECH  COUNCIL  R.  Scholl.  Palo  Alto  CA 

ARVIDGRANT OLYMPIA.  WA 

ATLANTIC  RICHFIELD  CO.  DALLAS.  TX  (SMITH) 

AUSTRALIA  Dept.  PW(A.  Hicks).  Melbourne 
BECHTEL  CORP.  SAN  FRANCISCO.  CA  (PHELPS) 

BELGIUM  HAECON,  N.V.,  Gent 

BETHLEHEM  STEELCO.  BETHLEHEM.  PA(STEELE) 

CANADA  Mem  Univ  Newfoundland  (Chari).  St  Johns:  Surveyor.  Nenninger  & Cheneveri  Inc..  Montreal 
CF  BRAUN  CO  Du  Bouchet.  Murray  Hill.  NJ 

CHEVRON  OIL  FIELD  RESEARCH  CO.  LA  HABRA. CA(BROOKS) 

CONCRETE  TECHNOLOGY  CORP.  TACOMA.  WA  (ANDERSON) 

CONRAD  ASSOC.  VanNuysCA(A.  Luisoni) 

DRAVO  CORP  Pittsburgh  PA  (Giannino):  Pittsburgh  PA  (Wright) 

NORWAY  DET  NORSKE  VERITAS  (Library).  Oslo 
EVALUATION  ASSOC.  INC  KING  OF  PRUSSIA.  PA  (FEDELE) 

FORD.  BACON  & DAVIS.  INC.  New  York  (Library) 

FRANCE  Dr.  Dutertre,  Boulogne:  L.  Pliskin,  Paris:  P.  Jensen.  Boulogne 
GLIDDEN  CO.  STRONGSVILLE,  OH  (RSCH  LIB) 

GLOBAL  MARINE  DEVELOPMENT  NEWPORT  BEACH.  CA  (HOLLETT) 

GRUMMAN  AEROSPACE  CORP.  Bethpage  NY  (Tech.  Info.  Ctr) 

HONEYWELL.  INC.  Minneapolis  MN  (Residential  Engr  Lib.) 

HUGHES  AIRCRAFT  Culver  City  CA(Tech.  Doc.  Ctr) 

ITALY  M.  Caironi.  Milan:  Sergio  Tattoni  Milano:  Torino  (F.  Levi) 

MAKAI  OCEAN  ENGRNG  INC.  Kailua.  HI 
JAMES  CO.  R.  Girdley,  Orlando  FL 

LOCKHEED  MISSILES  & SPACE  CO.  INC.  Mgr  Naval  Arch  & Mar  Eng  Sunnyvale.  CA;  Sunnyvale  CA 
(Rynewicz);  Sunnyvale.  CA  (Phillips) 

MARATHON  OIL  CO  Houston  TX  (C.  Seay) 

MCDONNEL  AIRCRAFT  CO.  Dept  501  (R.H.  Fayman).  St  Louis  MO 
MEXICO  R.  Cardenas 

MOBIL  PIPE  LINE  CO.  DALLAS.  TX  MGR  OF  ENGR  (NOACK) 

MUESER,  RUTLEDGE,  WENTWORTH  AND  JOHNSTON  NEW  YORK  (RICHARDS) 

NEW  ZEALAND  New  Zealand  Concrete  Research  Assoc.  (Librarian).  Porirua 
NEWPORT  NEWS  SHIPBLDG  & DRYDOCK  CO.  Newport  News  VA  (Tech.  Lib.) 

NORWAY  DET  NORSKE  VERITAS  (Rorent  Osio;  I.  Foss.  Oslo:  J.  Creed.  Ski:  Norwegian  Tech  Univ  (Brandtzaeg). 
Trondheim 

OFFSHORE  DEVELOPMENT  ENG.  INC.  BERKELEY,  CA 
PACIFIC  MARINE  TECHNOLOGY  LONG  BEACH,  CA  (WAGNER) 

PORTLAND  CEMENT  ASSOC.  SKOKIE.  IL  (CORELY):  Skokie  IL  (Rsch  & Dev  Lab.  Lib.) 

PRF.SCON  CORPTOWSON.  MD  (KELLER) 

RAND  CORP.  Santa  Monica  CA  (A.  Laupa) 

RAYMOND  INTERNATIONAL  INC.  E Colle  Soil  Tech  Dept.  Pennsauken.  NJ 
RIVERSIDE  CEMENT  CO  Riverside  CA  (W.  Smith) 

SANDIA  LABORATORIES  Library  Div.,  Livermore  CA 
SCHUPACK  ASSOC  SO.  NORWALK,  CT(SCHUPACK) 

SEATECH  CORP.  MIAMI.  FL(PERONI) 

SHELL  DEVELOPMENT  CO.  Houston  TX  (E.  Doyle) 

SHELL  OIL  CO.  HOUSTON.  TX  (MARSHALL) 

SOUTH  AMERICA  N.  Nouel,  Valencia,  Venezuela 
SWEDEN  GeoTech  Inst:  VBB  (Library),  Stockholm 

TRW  SYSTEMS  CLEVELAND,  OH  (ENG.  LIB.);  REDONDO  BEACH.  CA  (DAI) 

UNITED  KINGDOM  Cement  & Concrete  Assoc  (G.  Somerville)  Wexham  Springs.  Slou;  Cement  & Concrete  Assoc. 
(Lit.  Ex),  Bucks:  D.  Lee,  London;  D.  New.G.  Maunsell  & Partners,  London:  Library,  Bristol:  Shaw  & Hatton (F. 
Hansen),  London;  Taylor.  Woodrow  Constr(0!4P),  Southall,  Middlesex:  Univ.  of  Bristol  (R.  Morgan),  Bristol 
WATT  BRIAN  ASSOC  INC.  Houston.  TX 

WESTINGHOUSE  ELECTRIC  CORP.  Annapolis  MD  (Oceanic  Div  Lib.  Bryan);  Library.  Pittsburgh  PA 


54 


WISS.  JANNEY.  ELSTNER.  & ASSOC  Northbrook.  IL(J.  Hanson) 

WOODW  ARD-CLYDE  CONSULTANTS  (A.  Harrigan)  San  Francisco 

AL  SMOOTS  Los  Angeles.  CA 

BROWN.  ROBERT  University.  AL 

BULLOCK  La  Canada 

F.HEUZE  Boulder  CO 

CAPT  MURPHY  Sunnyvale.  CA 

GREG  PAGE  EUGENE.  OR 

R.F.  BESIER  Old  Saybrook  CT 
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